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1 . Introduction 


Many nonlinear dvnamical systems involving density waves 
can be represented by a system of partial differential 
equations: a nonlinear equation of wave propagation and an 
equation of evolution for a certain dvnamical variable. The 
dynamical variable may be a field as in the Schrodinger 
equation for solitons, or mav be a fluid velocity as in the 
Navier-Stokes equation of motion for compressible turbulence. 
The general study of coupling between the kinetic energy of 
turbulence and the wave energy belongs to the topic of acoustic 
turbulence, * »— , while the escitations of turbulence by wave 
motions cause a density-induced turbulence. The latter is 
treated here, because it occurs more frequently in the 
atmosphere in view of the general presence of the internal 
gravity wave and the mean density gradient. 

A kinetic method is developed in order to include those 
collective phenomena which are otherwise not evident. By group- 
scaling, the master equation is decomposed into a hierarchy of 
scaled equations for the velocity distributions of decreasing 
order of correlation times, representative of the macro- 
evolution, the transport property and the relasation. A cluster 
of distributions participates into the role of relasation for 
the approach of the transport property to equilibrium. I f the 
cluster loses the individual identity of velocity distributions 
bs' moving together and behaving like a clump, an effective 
medium can be assumed and a closure is found. 


j 
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By means of propagators that are modified to include the 
damping from the effective medium, the eddy collision is 
calculated as an operator, and the kinetic equation for the 
macro-distribution is derived. It contains a memory 
representative of the collective behavior. 

Based upon the kinetic equation. the transport theories 
for the eddy diffusivity and eddy viscosity are developed. Thev 
calculate the transport functions to govern the following 
processes. A mean density gradient or a large scale gravity 
wave may escite all the nonlinear mechanisms that control the 
density waves, i.e. the coupling mechanism for the excitation 
of turbulence and the cascade transfer from large toward small 
eddies. By the more efficient coupling mechanism, the density 
waves feed energy into turbulence that cascades down the 
velocity spectrum. 

The transfer function is found to have two memory-loss 
functions: one from the eddy collision as an operator, and the 

other from the Lagrangian correlation between the velocity 
distribution and the fluid velocity. The two memory-losses 
formulate a direct cascade toward small eddies and a reverse 
cascade toward large eddies. respectively. The former is 
relevant to the coupling subrange. 

The spectral intensity in the coupling subrange are found 
to have a k“- law for the velocity and field (i.e. pressure 
gradient) fluctuations. and a law for the density 

fluctuations . 

The formation of sawtooths^ in compressible turbulence can 
i 
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be seen in numerical computations . yielding a 

intensity k~“ , in deviation from the Kolmogoroff 
incompressible turbulence. 


spectral 

1 a mP ' ® o f 


2. Basic equations for the description of the microdynamical 
state of density-excited turbulence. 

The hydrodynamic variables are the fluid velocity u. the 

A y\ 

density and the pressure p . We write the density 


N < t , X > 


f < t , X > / P 

- CO 


< 1 > 


that is normalized by means 
the thermal speed 


of a constant density ? 

<x> 


and write 


A 

C 



(2) 


A fluctuating variable is denoted by the symbol < >, with 


( > 


/V 

+- < > 


< 3 1 


It consists of the superposition of the ensemble average 
' '*<> and the deviation from the ensemble average, called the 

fluctuation < > . For such a decomposition, we may use the 

operators A, A, such that 

n — yv 

A d A -h A , < u > 


4 
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where 3=1 is the unit operator. The decomposition of a 
fluctuating function into an average and a fluctuation. as by 
is known as the Reynolds decomposition in turbulent flows. 
For a proper selection, we first examine the following models 
for the description of the microdynamical state of turbulence. 

( a > full Navier-St okes model 

The Navieon-Stokes system consists of the equation of momentum 

d, N + V. /V 14,. OL. — f\/ E . ( 5 ^ 

t i / / * 0 

and the equation of continuity 

\ N + V ft tJ-. - o . ( 6 ) 

t J </ 

n 

The driving field E is due to the pressure gradient and is 
given by the relation 

E = - N 

- “* ( 7 ) 

A A 

between E and N. The differentials are 

'b/K , V s *Ax • 

and like other operators apply to the functions which follow. 

( b > acoustic turbulence model 

For describing the interaction between turbulence and 
finite amplitude density waves, we can transform the full 
Nav i e r ■■ St okes system in model <a> into the following simpler 
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momentum equation 


(\ ' *■*)* = i 


( 8 ) 


and the equation of propagation of density waves 


(>£ Z ~ v 7 * c z v)N = st , 


< 9) 


with 




-'I 5 


V7-- W «•«. 


< 10 ) 


/» /» 

The relation between E and N remains to be (7). 

w* 

The equation of propagation is obtained by cross 

differentiations of (5) and (6) with respect to t and x and by 

«■* 

a subtraction 

It is to be recalled that the simple form of momentum 

A* 

equation (8) is obtained by an elimination of N between the two 
equations (5) and (6) of the full Navier-St okes system. The 
same simple form of momentum equation is valid for 

incompressible turbulence as well with the divergence-free 
condit ion 


V . LC ^ o . 


( 1 1 > 


j_cj_ Ri emann ' s model 

In order to show the coupling between momentum and density 


) 
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in a concrete manner, we transform the Navier Stokes system 

7 

into the Riemann system in the form : 


(\+- CL-vJu - - c vw <i2> 

('b +U,-v)w*.-CV'U. . <13) 

La. ' M. M 


TheRiemann variable 


^/ — C ^*L A/ (14) 


relates w and N, with a constant speed of sound c in isothermal 

8 

gas and a variable speed in adiabatic gas. 

Upon multiplying ( 12' ) and < 1.3 > by u and w, respectively, 
and averaging, we find the equations of evolution 


a (^-h u,- v) - w <i5> 

4 -h &.v)w)* -W, (16> 

for the kinetic energy and the potential energy in homogenous 
turbulence. The equations are coupled by the function 


W s <171 

with a net zero balance. 

With the supply 'by external source, the wave energy feeds 
an amount W to the kinetic energy. The latter then undergoes an 
internal cascade transfer from large into smaller eddies, to be 


i 



ultimately' dissipated by viscosity. 

(cT> kinetic model 

It is not difficult to verify that the two equat 

and <61 of the Navier-St okes system are the first two 

9 

of the master equation, written as 

(V ' l )\( t 'i' v ) = 0 

for the distribution function in the form 

Nfa*) $L V ~ 

as the condition of the normalization, i.e. 

- A 

The differential operator in the phase space t.x.v is: 

«%• 

L(w; s v-y + eft. ,*;•) 

with d - 'V ^ s’ . 


<-ed wavfi - k i ne t i c model 

The present model uses the master equation ( IS 5 


ions (5) 
moments 

< IS) 


< 19 ) 


( 20 ) 


<21 ) 


with the 


condition 
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to replace the momentum equation (8> of model while 

retaining the equation of propagation ( 9 > and the relation < 7 > 

for defining E-field. 

•v 


Being homogenous and having the velocity v as an 

independent phase variable, the master equation has lesser 

nonlinearity. It is most equipped to treat the collisionless 

transport properties for the wave -induced turbulence. It 

provides a mechanism of collisionless damping as the result of 

the interaction between the velocity distribution and waves. 

This property cannot be derived from the hydrodynamic approach. 

Among the kinetic methods, the wave-kinetic model is preferred, 

because it describes more explicitly the interaction between 

finite amplitude density waves and turbulence. 

The transformation of the Navier-St okes equation into the 

master equation and the subsequent derivation of the kinetic 

equation of turbulence have been attempted by Mon in 10 and 

Lundgren . 1 1 But thes encountered the difficulty of kinetic 

1 2 

hierarchy like in Bogoliubov's theory. To break the hierarchy, 
we introduce a method of decomposition into many scales. 


3 . Multi-scale group decomposition 

The hydrodynamic equations and the master equation 
discussed above describe the microdynamical state of 
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turbulence. They contain too many minute details which are not 
suitable and necessary for the statistical treatment. We 
introduce a course-graining procedure called group-scaling by 
writing 

rv ° 

A - A - A' (- 2 . 3 ^ 

A ' = A 1 * A” , < 23b ) 

as an extension of the Reynolds decomposition < 4 ) . The three 
groups 

A° , A ’ , A" ( 24a ) 


are called the macro-group, the micro-group and the submicro- 


group, 


respectively. while A x is called the submacro-group. 


The groups (24a) and the following groups 

, .A 1 , A W , . . . 

must conform to the scaling differential 


(24b) 





(25) 


of decreasing duration of correlations, indicating a decreasing 
coherence . 

The three groups (24a) represent the three processes of 
spectral evolution, transport property and relasation. The 
latter governs the approach of the transport property to 
equilibrium. The decomposition gives the possibility of 
exploiting the property-' of local quasi-homogeneity between any 
two groups. This assumption implies that the interaction among 
the groups are restricted to nearest neighbor group-pairs. 

The fluctuating groups, as formed by the scaling operators 


j 
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(2 i), are denoted bv the superscripts 


o i (z ) 

(>,(),(>.... ( 26 a > 


and 

o 

(>,(>•,(>" ( 26 b) 

in the open sequence and the closed sequence, respec t i ve l S' . 

For completeness we should add the cumulative operators 


A 

o 




A’ 


A’ 


giving the cumulative groups 

( > , ( > , ( > 

<» I Z 


A 7 


( 27a) 
( 27b ) 
(27c) 

( 28 ) 


The deterministic transport properties and transport 
functions, as shaped by these fluctuating groups, are denoted 
by the superscripts 




( > 


< ) 


Cl? 


< > 


C'l 


( ) 


V 


( 29 ) 


with square brackets. 

The fluctuating groups may be Fourier transformed to find 
their cumulative spectral distributions, as follows 


JL 

Z 



( 30a) 
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i(E' Z } - fjLk" F(k") 

A 

= J°kk~F(k"'), 

k" £ 

with the inequalities 

k > k " > k ^ k' 

among the wavenumbers. They are independent variables or 
variables of integrations. The limits of integrations are the 
demarcation wavenumbers of the three groups. The cumulative 
spectral distribution obviously yields the spectral density 
F ( k > of field fluctuations by a simple differentiation with 
respect to k . 


< 30b) 


<300 


4. Kinetic treatment of turbulence 

4 . 1 Effective rtFdtrnn approximation for the closure of k inetic 
hierarchy 

In the statistical mechanics of many particles, the 
master equation generates a hierarchy of man>— part icle 

distribution functions. Bogoliubov closed it at the triplet 
distribution in an arbitrary manner. ” Analogously, our master 
equation is expected to generate a hierarchy of scaled 

J 
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equations by the group decomposition. We may choose a 
closed sequence ( 24a > or an open sequence (24b'). The latter 
contains more information and is investigated here. 

The scaling of our master equation by the operators of 
open sequence yields the scaled equations 


L ° 1 { * = - L * f 4 z ° - z 

(k+Ljf ,-L i (f\f)+ z 1 - z;... 


for the distribution of open groups, with the differential 
operators 

L . L , . . . (34) 

o 


from (21). The scaled equations have their eddy collisions 

Z'-AL'f - 

z‘*-alY~ 

Z WlT- 


with the collision operators 



(36) 


Through collisions each group i s coupled to the next one 
by developing a hierarchy. Physically speaking^ (32) determines 


J 
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the evolution of the velocity u and of its kinetic energy, and 
< 33 > determines the transport property by <.35b’) . Finally, the 
cluster 

f" = /** * . . - < 37 ) 

organizes a relaxation process for the transport property to 
approach its equilibrium and is called the relaxational 
cluster . 

Two methods are available for the closure. Firstly, the 
equation 


L 7 * 

without collision may serve as a closing equation. It will 
calculate the collision by writing 

Z l a Z' < 39 > 

and 

z'--A’Ly. 

The early cutoff neglects the collective behavior of 
long range. Secondly, we realize that the problem of hierarchy 
can be found in many nonlinear dynamical systems. The familiar 
ones are the hydrodynamics of long chain molecules and the 
suspension of particles in Brownian movements, ^ For these 
systems, the concept of an effective medium*^ that offers a 


i 
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Darcy damping has been introduced. Likewise, Heisenberg 
postulated an effective viscosity in his theory of 
turbulence. Hence as the second method of closure, we shall 
assume that the relaxational cluster of many distributions 
represents in its ensemble a porous medium that is frictional 
and presents a similar Darcy damping. We shall determine this 
damping by a self-consistent theory of collective phenomena. 

To elucidate this physical picture analytically, we take 
the moment of (33) to get the momentum equation in the form: 


i ri 


(\ “r v ) a± * - ^ + E + j 


The hydrodynamic friction is found as 


t ■ 1 


A v v ~2L 


= V 

i. 

= - y 6t / 


(.41) 


(42) 


by (391 and (40). It is proportional to the coefficient of 
damping g . Although it assumes the same role as the Darcy 
damping, it finds a self-consistent kinetic basis in (42), 

For the role of relasation, it is more important to 
retain the collective behavi' or than to discriminate the phase 
individuality among the distributions in the cluster. Hence we 
can replace the cluster by an effective medium with the property 
(42), and write the approximation 


i 



This simplifies (33) into the following 


(\ +£ ±)f L --^G*P.- z °. 

with the differential operator 

£ ■ L i + ^" 7 < 45) 

5 . 2 Local q uasi -homogeneity 

The scaling differential (25') permits the assumption of 
quasi-homogeneity among groups. Thus a low-order group of long 
correlation time can be assumed to -vary more slowly than does a 
higher-order group of shorter correlation time. BY applying the 
property of quasi-homogeneity at various levels, we find that 
the interactions between the nearest-neighbor groups are the 
most pronounced. Under this circumstance, we can simplify (32) 
and (44) into 


(yUf = - - z m 


5 . 3 Macro-kinetic equation 


The integration of (46b) will be made by the intermediary of the propagator 
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satisfies the differential equation 

(\^)U(t,tJ =o. 

The omission of Z in <46b1 implies 9 


1 -*) . 


By assumption (48), we may use the average propagator to 
calculate 


f i =- fW A L {( , 


(49) 


.1 1 

with^,=L ■ b Y a substitution into (35b>, we calculate the eddy 
collision 


z”= - A L i j i 

= ep- 7 {f} 

s <50> 

with the collision coefficient 

f/~\ r ~t -*co 

£> = J <h=(L'd*, v jl!!(t,t-ijL(t--cjy 

= JD 1 ' 1 - 5 <51' 


and the eddy diffusivity 
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T$ 1 = J Vi e (t-*)) 


( 52) 


By the scaling differential 


07 


O T 

the upper limit t has been put to o o without altering the value 
of the integral. The asymptot i c i tv renders the diffusivity 
deterministic . 

Finally upon substituting (50) into <46a>, we derive the 
kinetic equation for the macro-distribution in the form 


h+ULy.-L-f + f 


( 53a) 


that can be rewritten as 


C\* L + L°)f=-L‘f + V -fhn 


( 53b> 


by (51). The macro-kinetic equation possesses a memory as 
belonging to the non-Markovi an behavior through the operator 


> C ^} 


5. Derivation of the hydrodynamic equations of turbulence from 
the kinetic equation 

By taking the first moment of the macro-kinetic equation 
(53) and with the condition (22), we obtain the momentum 


7 
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equat ion 


t * * vu. 


= £% J 


The hydrodynamic friction takes the form 


j”=tw ' 

or 

r-jh* ?f} - 


< 54 > 


< 55a > 


( 55b ) 


by (51>, 

Since the macro-kinetic equation (53') is an equation 

espl i it in f° without the involvement of f ’ , its moments will 

o 

generate u and not u' , This ms a consequence of the irrever- 
sibility of the macro-kinetic equation. This property has been 
observed in deriving the momentum equation (54). 

The equation of kinetic energv is obtained upon multi - 

o 

plying (54) by u and taking an ensemble average. It has the 
form 


Co3 r»7 Co 7 

s = w - t , 

with the following transport functions: 

( i ) The function 

is the rate of change of the kinetic energy. 

( i i ) The transfer function 


(56) 


( 57 ) 


1 
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- i " C P “] / O 0 \ 

T = " / U, • J ) <5S> 

governs the cascade transfer across the spectrum. The direct 
cascade is a transfer toward the large wavenumbers, and the 
reverse cascade is a transfer toward the small wavenumbers. 

( i i i ) The coupling function 

t/ 1 . <u*-e*> 

- < 59 > 

governs the coupling between the kinetic energy spectrum and 
the density spectrum. It represents an excitation of the kine- 
tic energy of turbulence by the density rarefaction, i.e. 


fo7 

W?o. 


6. Transport theory of diffusion 

6 . 1 Laarannian correlation 
The diffusivity dCO 

that is basic to the macro-kinetic 
equation is defined by (52) as the time integration of the 
Lagrangian correlation 

<e ftp Tifty-ft E'Ct-vi) 

of E' -field fluctuations, where '1$ is the propagator for the 
evolution of E'(t-T') along the trajectory that is perturbed by 
the differential operator (45"). 

By Fourier decomposition of E'(t,x^ into the frequency " 
and wavenumber k", we have 

L -k"*J 


*** JJ ' V*. 


( 60 > 
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The Lagrangian correlation has the form 

(ejt c ) e(t-z ) > = J'fak'). 

s uk -Jlfc) \ 

* * * (e ~ ~ / ' 




(61 ) 


where yC<T ) is the path-length in the time interval T . The 
damping yM 

as a differential operatorhas its Fourier f o r 

O'] 

Y . The first factor 


m 


i'W; * / (e'Hki k."iy 


(62) 


in the right hand side is the spectral function for the field 
fluctuations, and 



(63) 


is the factor of truncation in a Fourier decomposition that is 
truncated within a time interval 2T and a length interval 2X in 
three dimensions . The spectral function is followed by a number 
of esponential functions, called the components of the orbit 
function . 

The component 

_ i(o3 k'vjz 

S? -e. ~ (64) 



represents the linear streaming at velocity v in an unperturbed 
trajectory. It defines the Eulerian correlation 


a 
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(£^vj£fH^vrj) = JJjco''ctk" Jty'k •'/ UuX'zv ) , 


(65) 


and forms the diffusivity 

PI rr ... i'l 


VkVfjz U» 7 favJ 
= f p'j-W'J'^k-) r("«"-Av; 


in weak turbulence . 
The component 


UKz) * ‘ 


represents the effective medium. By 

< - W', 

•V X 

the esponential function 




•e ~ _ y = ^ 


can be written as the product of two components 

l(fc) a <«. 

(ii — 


( 66 ) 


<67> 


( 68 ) 


(69) 


(70) 


(71) 


They represent the nonlinear streaming and the diffusive 
relaxation, respectively, along a trajectory that is perturbed 
by turbulence . 

We are dealing witn strong turbulence, so that the linear 


J 
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effect by the orbital component 



for streaming is negligible, escept when calculating the 
collision < 51 > . 

By collecting the components into the orbit function 

Uk w ) -uk^j i (k'j-cj < 73 > 

• Cl]~ - ' 

we can write the Lagrangian correlation in the form: 

(tfcj *) lift, (r-K ) t'ft-tj) = jjfl 4 \k") 4v(k l t) < 74 > 

by ( 6 1 > , with 

* %{E{k"j B'(.i ?> 

w w W ^ v "> 

and < 75 > 

^ £ (?/a/ ■ 

The diffusivity < 52 > , rewritten as 



has a correlation time 

*** *0 V% w / ) 

Here and in the following the integrals without limits are 
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understood to estend from -co to®. 

6 . 3 Path dispersion 

The evolution of the orbit is governed by the dynamical 
equations in the form 


i *~ ( V- u(t ) , Ut } , 

- 1 4t. ~ 1 


( 78 > 


with initial conditions 


ic(t = t ) - X , = V / I fr - E it, x) 


( 79 > 


or, in the alternative form 






< 80 ) 


with the conditions 


^ (t-o) ~ 0 , \/fc = <Dj ^ 0 


at "C = 0 . 

The path 


<81 > 




< 82) 


can be decomposed into an average 

ji-0) = fE -C Z 

«V 

and a fluctuation 


( 83 ) 


s 
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< 84 ) 


In stationary turbulence, the t-dependence can be neglected. By 
(8D, the path is independent of v. 

From (84> we calculate the variance 

(hjl( t)) =2 / T ‘"i - TV +t T 'y(e/V < 85 > 


the diffusivity by path dispersion 

- -C J X 4t'(x--z') (e{t) e(^)) < 
and the integral 

i 4 r * (5A ’) s y 


^ (t T,_ 2 T T + 


The asymptotic values are: 
( i > for large X 

<A>) AA = T x * 




K (t] ■ T>(x-Kb) x 2 

=• 

i J> 

with 

a4.feV.-3> ; 

3> 3 w ^ * 


( 86 ) 


(87) 


<88> 


(89) 


<90> 


(91 > 


l 
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< i i > for small X" 


K/-^ ? t ^ (z&iGj) 

VM 

K 0 (t") £ 7w^ 

0 x /f 


( 92 ) 


< 93 ) 


< 94 ) 


with 


'Jw 


* 3 JL 


i 




/>] 


( 95 ) 


6 . 3 Orbital component 

Introduce the probability of retrograde transition 

fo,P 

w / 

with the condition of normalization 


< 96 ) 


M fop = i 


< 97 ) 


to calculate 


A"1 / f i a Ck '■£ , n) 

l (fo)-fo * “ >=}*{ + 


< 98 ) 


The dispersion of path fluctuations is governed by the 


following equation of transition 




( 99 ) 


where is the diffusivity for the dispersion of path 
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rewritten in the form 


l"l 


K ( (x) = T JV v '( t- t) , 


( 100 ) 


from < 86 > . 


By -Fourier decomposition with the formula 

bfcjl = -! — M 

1 " (lx)' 1 ~ 1 - 

we transform <99> into 


< 101 ') 


OT 


\ £7 = - k 'V •• K. (?) fa k "} / 


< 102) 


and obtain the solution 


f> h k "1 = -i— *4 f- n ph’ kX o 

~ (totf IL l 


( 103) 


The coefficient is f° un d by the condition ( 97 ) , when we put 

k” =0 in <101>. 

«V 

By comparing <98> with <101> and using < 103 > we find 




■ r- iv-'i? * x x 


- % 3 t } ) , 


■0 ~J 

R3 j 


< 104 > 


with 


^o L = _L le A : D £ -u fe X) 

D 3 - - - 3 * < 105) 

The asymptotic value <90> at large ~C is taken, in view of 


i 
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the scaling differential 

T > X 
c 

The diffusivity dR as belonging to the highest group, 
regulates the relasation, and will be cal led the relasational 
diffusivity . 





6 , 4 Orbita l c omponent 

BY the same procedures as for (104), we calculate 


V- ;t)= Lk ~~ h) 
= “K- 


i \ *** ' 


< 106) 


wi th 


So = f " ■■ 




( 107) 


The asymptotic value (94> at small X is taken in view of the 
scaling differential 


T<I 


6 . 5 Collision diffusivity D^ 1 ^ 

By collecting all the components 



ik'-v x 
-e 


V w, 

E 

Y w 

- o b -t 
e 

< 108a) 

ixY^yk") ~ 

- % 3 ~c 3 

-e 

< 108b) 
( 108c) 


3 
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< 108d> 


of the orbit function 


i(r f k'y) 


from <73’), < 72 > , (lOt'), and <106>, 


we determine the correlation time 

r c \ s V J 


w 


= J 4> v (Kt>, v) 'Ly(k'P) CA/ V 


< 109) 


for the collisional diffusivity, which can now be written in 
the form 


• and 

k 

a tensor and a trace, respectively. 


( 1 lOa') 


( 1 10b) 


6 . 6 Relaxations! diffusivity D ^ ^ 

From our basic kinetic ' system ( 46 > for fo and f£, 

that was used to derive the kinetic equation (53~> and the 
diffusivity D it is seen that the submicro-group f" has no 
esplicit governing equation like <38>, but is represented by an 
effective medium without identification of individual velocity 
distributions. This means that the collisional diffusivity 

and the relasational diffusivity of groups 

M 1*1 

D . D , . < 111 > 

have to be treated differently. While the collisional 


j 
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di ffusivity is calculated by including all orbital 

- a * 

components, as being based upon the propagator, the 

relasational diffusion tv of various groups do not have their 
individual kinetics and propagators, because they all form 
part of the same effective medium. For this reason, the 
dif fusivities ( ill ) are self-generating in the symbolic form 


M t'A. 07, 

D = D JD ■ } 




bv selecting the diffusive component * o 

* 

function . 

More esplicitly, (112) is 

~ 07. 


( 112 ) 


f the orbit 


with 


m tf w cr) , 

k" 


(113) 




M -i 


( 114 ' ) 


Upon substitution, we obtain 

4' 1 ' 1 

k' 


-1 


= fr(f)f^'F(rJ (t k"' j> M )~ . 

The use of (110), (77) and (.1080 has been made 

r /"i . J>7 


( 115) 


The d i f f us i v i t i es 


and D 


have an identical 


integral structure, but with contingent limits of integration 


in 


( k" , oo > and ( k , oo ) 


i 
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r^j r*>] 

respectively, so that we can identifv Er as a diffusivity D 
in the Integra^ making (115') an integral equation. For its 
solution, we differentiate (115> with respect to the lower 
limit of integration k", as denoted by an upper dot, and obtain 
the differential equation 

i> M = - f raj Fjk ") /(j k" 


This is solved to give the diffusivity for relaxation in the 
form 



( 117 ) 


M 

The relasational diffusivity D forms the last link of the 
chain of many d i f f usi t i v i t i es . The approximate formula 

provides with the necessarv relaxation process for the approach 

[/] 

of D to equilibrium and for the closure of the sequence of 
d i f fusivi ties, 


7 . Enhancement of turbulence by density waves 

The momentum equation (54) can be written in the following 
alternative form: 



( 118 ) 


with the differential operator 

.C'l 


£ s ^ v + Y 


( 119 ) 
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o o 

The nonlinear term represents the steepening of sawtooth 

waves. The steepening cannot increase indefinitely to give rise 
to a discontinuity in view of the damping by v . This damping 
has a kinetic origin as due to the interaction between the 
waves and the distribution function. 

By considering the left hand side of <118> as the total 
time derivative, we can integrate to obtain the macro-velocity 


it° =■ [ J* E°( i--*) , 

<0 


< 120 ) 


and the diffusivity 




( 121) 


S') 


The asymptot icity is obtained on the same basis as for D . The 
evolution operator”^ determines the Lagrangian evolution of u 4 
under the differential operator (119). The same earlier 

argument has been used for approximating T$ bv V‘ . The operator 

_ u w 

~\Jt , being of hydrodynam i c origin, does not involve v in the 
H CoJ 

perturbed trajectory, so that D has the same structure as 

f/] co . 

D from <H0a), but D is independent of v and has different 
limits of integration. Thus we can write 


1 > H = F/k'/ \\k') , 

V A 


< 122 ) 


with the correlation time 


J 
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-rjl') = p T l(k' /X .^ 0 ] . 


< 123) 


The orbit function 


^(K^^ro) = ^y(kp) Vhty ^\k'-c) 


and its components 


^y(k'-u) S 

4/^ s 

^ (K*) s e 


-yf'7 

-/ T 


— T ^ 
G»1 ^ 


, . Mi J 
T 


< 124) 


< 125a) 


( 125b) 


( 125c) 


are modified forms of <73) and (108). Here the diffusion time 

r P~ l 

is H) , such that 


^5 J> 


Mi 


r 

>=0 


and the diffusivity 


3> 


M 


v-0 

has a correlation time 


= /4" fVi W( , 


p\A = / ^ 


v -0 v o 

with the orbit function 


V=0 




M 


< 126) 


( 127) 


>=o y C13 

The use of (105>, <110a), <109) and <73) has been made. 


( 128) 


< 129) 


i 



By comparing the trace of the diffusivity tensor < 121 > 


with the coupling function (59) we find 


W = 3 D - 


< 130 > 


with 


tP 1 a tUoi 


< 131 > 


This diffusivity governs the enhancement of turbulence from the 

o o 

coupling between E and u . In the spectral sense there is an 
enhancement of the velocity spectrum at the expense of the 
density spectrum by an araout 


8, Transport theory of cascade 

8.1 Two memories of the transfer function 

In sections 4 and 5 , we have derived the macro-kinetic 
equation < 53 > and the hydrodynamic equation of momentum < 54 ^. 
The eddy collision is 



< 132 ) 


by (51 > and the hydrodynamic friction is 
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by < 55b ) . 

o 

Upon multiplying by u and averaging, we have the 
transfer function (58) which we rewrite as follows: 


-r c °1 / . ov r v r r®7 


< 134 ) 


where 



( 135) 


is a Lagrangian correlation. 

Like in <61 > we use the Fourier method to transform (135) 
and write 


(v X,\/) - Jjjt' A'Ct* ) *f (k ) 


(. 136a) 


with the Fourier component 

i.\') - 

By taking the moment, we have the spectrum 

v i ^\k')= t\k 1 

L -V* 

and the energy 



with the definitions 


( 136b) 


( 137a) 


< 137b) 


4 
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't ft *} = /<#' - <“•”!> . 

Bv substituting for the Fourier form of 

2> C ' 7 and f C ° ] 

from (110a) and <136a), we transform (134) into 

-f°\ - ffol'i** J-'J (kip r N*?(k:k'J 
1 JJ — ~ 'b ~ - 


< 133a) 


( 138b) 


( 139) 


For the sake of abbreviation, we introduce the partition 
function 




< 140) 


The dependence of 


o pI 


W). W), % (V 


< 14la) 


on 


,oi _.M m w ro 

t y , y , y / 


( 141b) 


is understood by their definitions. 

In order to simplify the expression (140), we consider the 
following properties: 

( i ) In strong turbulence, the streaming effect bv v in h 

-w. V 

is negligible, escept with the derivatives 


a 
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YU<w -- - *; K i,m I 


VcO 


= -J- k' z s. 

5 J* 


( 142a) 


and 


< 142b) 


t'3 


\\ - - f k' (f- * 

and with similar derivatives of h<k"). 

M 

< i i Since T is finite and real, and so is 

the integrand in <139"> must be even. This rules out the 
contribution by ll. 

< i i i ) With the partial integration with respect to v, we 

have 






( 143) 


(iv) The transfer function controls the transfer of energy 
across the spectrum. Consequently, the integrand of <139> must 

retain v . , 

4 

By means of these properties, we can simplify the 
partition function <140>. We write 


1 V- v V y * v J*-'! 

and calculate the two components 


and 


(rfprf** 


< 144) 


( 145a) 


145b) 
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The two memory-loss functions are: 


L / /L \ - </ V-0 - V=<3 


mmj -= u <w| 5* • 

^ U </ A - - v -o 3 - * vzo 


and 


( 146a) 


( 146b) 


The loss of v in the orbit function makes the memory-loss 

v% 

functions only dependent on the wavenumbers k', k" as scalars. 
It is seen that < 144 ) describes the partition of the macro- 

yf el 

energy of turbulence G <k') in the k'-space, while accounting 
for the memorv losses along the two orbits, of orbit functions 

ift'j and i(n . 

x. v_ 

By the two memory losses, the transfer function (139) can also 


be divided into two parts, as follows: 

roi r»7 

T = <T ) 

<*, 


M 

<T X. 


with 


2 V :k \ 


( 147) 


< 148a) 


( 14Sb) 


M. 


by (145). 

The transfer function i s energized by the macro-energy (k) 
and is diffused by the micro-field fluctuations of intensity 
sMk"h Since the process takes place in the Lagrangian space, 
the two orbit functions cause two memory losses. 


a 
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8.2 Direct and reverse cascades 

By noting the scaling differential (.25) and by making use 
of the properties of the partition function mentioned above, we 
reduce (146a) and (146b) into 

[HjWJl = - vjo) ^ 

^ v,oJ ^ . 

We can see that the loss of memory is effectively caused 
by the orbit function 


V- oj 


( 150a) 


as the result of the perturbations of the trajectory carrying 
the field 

Xlfat-T) e'(i- r) 

c.3 

for shaping the diffusivity D . The other orbit function 




< 150b) 


merely serves to provide a macro-gradient of a spectrum 

fc" Z%') . 

by the transformation 


i> l 'L(k r ) 


= - k' l -C Z v=oj • 

V= o ** 

to* 


( 151) 


The substitution of ( 149 > into (148) leads to the transfer 
function in two components: 


i 
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(■f') -.jjjLk'dk-J'frn' & (k-)jHxU(k:z,,-oj 

(t*) -jiuvi* 'f'(k')jy* t *( k ;^r=v 


with 


C' ] , f>J 

^ = LJX-9-C4 -O 


We write < 152) in the following concise form 


(t c °V - K" j k‘ 


I' 7 fel 


0 *i - ^ > , 

after a rearrangement by separating the variables k’ , 
governing functions are: the vorticitv function 

* M -<( v'iT? 


■ w *'(■*}> 

■ /•*■£' k' L f\k'j 


the eddy viscosity as a scalar trace 

0 


K 


■ 3 = JotV' ✓/ V/ 

’ f-pr Fjk'1 Cr%*l 


and the eddy damping 




1 PV/ */ r^fL»\ 


JO 


X'L 


*fl4k'k*Z(klCW. 

J& 


The modulation function 


( 152a) 
( 152b) 


( 153a) 

< 153b> 
k”. The 


<154) 

< 155) 


( 156 ) 


i 
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Q p \k") = Hr T 2 *[-%&]’ 

0 (157) 

£>? 

finds a modulation time W ■ 

The orbit function li(k",'C , v = 0^ represents a damping as 
the result of the perturbations of the trajectory. It 
determines the relaxation for the approach of the diffusivity 

C/1 

D to equilibrium. Likewise it determines, by (157), the time 
C'l 

7 -- of modulation for the approach of the eddy viscosity (155) 

to 

and the damping coefficient (156^ to equilibrium. 

The two parts (153a) and (153b) of the transfer function 
contain the second moments (154) and (156) of the spectral 
distributions with limits of integrations 

( o , k ) and ( k , ) , 

respectively. The two moments gain their importance for large 
and small values of k. respectively. Hence we can conclude that 
the two parts ( 153a > and (153b> of the transfer function 
control the direct cascade near a sink at large k and the 
reverse cascade near a source of small k. respectively. 

In his theory of incompressible turbulence, Heisenberg has 
proposed a transfer function in the form of a product of the 
vorticity function with the eddy viscosity5 His formalism was 
based upon the empirical concept of mixing-length. In addition, 
he proposed an empirical formula of eddy viscosity on the basis 
of dimensional consideration. Here we have developed a kinetic 
theory of transport for the transfer process. Our analytical 
formula of transfer (163a) by direct cascade supports his 
phenomenological formalism. Our theoretical result of eddy 

3 
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viscositv (155) is based upon the macro-kinetic equation and 
has undergone the process of relaxation, modulation and memory 
losses. It is obviously more complicated than Heisenberg's 
formula . 

En analogy with the derivation of the kinematic viscosity 
bs' the kinetic theory of gases, the derivation of the eddy 
viscosity needs the kinetic equation as an essential basis. 
This explains why a purely hydrodynamical method, e.g. the 
direct interaction approximation , 1 " 7 cannot produce an eddy 
viscosity . 


8.3 Modulation function 


.07 


The modulation function Q- (k") in (157) controls the 
eddy viscosity (155) and the eddy damping ( 156 ), and is 
calculated by the second moment of the orbit function of three 

components 

(ft) 

- ^{ •fa**)'} ,158 > 

of time scales as given bv the inverse of the frequency scales 



C'] 

dj "ViO. 

35 ' flj 


( 159) 


respectively. Each component could form its own individual 
modulation. An approximate evaluation of the second moment of 
(158) gives the composite modulation, as defined bv (157) with 


a time sea 1 e T 

M 


C'l 


such that 


3 
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T 


C'l -/ 








W2 
3 to 
3> 



< 160) 


is the resultant of the frequencv scales (159). The numerical 
coefficients are adjusted in such a way that the formula (160) 
gives the esact values of the three individual modulation 
functions at the three- scales (159). Thus the formula (160) for 
the composite modulation is an interpolation of three esact 
components . 


8 , 4 Identification ©£ effective medium 

f/J 

From (153a), we see that the eddy viscosity K describes 
the turbulent property of the medium of small eddies that carry 

the E' -field fluctuations. 

»%* 

On a phenomenological basis, we can conceive that this 
medium offers a friction 


7 = t\ v u. 

o 

to the evolution of u , and causes an energy dissipation 


( 161 ) 


7 V> - - K C %^T> 
= - K h 


( 162) 


to the energy balance in homogenous turbulence. Since the 
energy dissipation (162> is the result of interaction between 
the micro-group and the macro-group, we see here that the 
physical picture of cascade emerges. By pursuing the comparison 
of (161> with (42), we can identity 




} 


( 163a) 
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for the micro-group, and similarly for the submicro-group 




V 


z_ 


v". - k 1 jc M 


< 163b) 


Recall that the realization of a direct cascade requires a 
transfer of large wavenumbers in a transport of the gradient 
tspe. Both conditions are met by the effective medium that 
represents the cluster of high-order distributions. 

The identification of the effective friction eddy 
viscositv in (163) connects the concept of effective medium to 
the kinetic basis of cascade process in a self-consistent way. 


9 . Equation of state of turbulence 
9 . 1 Local nonlinearity 

The density-induced turbulence is governed by the 
following system of two differential equations: 

(\ * = ? 

(p£ -VC*v)ti - fl 

with the driving forces 


< 164) 


( 165) 


£ =■ - c 1 vJU, N 

" ~ <166) 

and 

/l - VV ■■ b) U-U. , 

< 167) 


from <S>, (9), (7> and (10), respec 1 1 ve l s - . 

system describes the m i c rods’nam i ca 1 state of the 


i 


The 
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interaction between turbulence and finite amplitude density 
waves. The Navier-Stokes equation ( 164) is obviously nonlinear 

A A 

by its advection u -Vu., and is also nonlinear by the driving 

A 

force E. The latter nonlinearity arises from the wave 
propagation by the definitions (166) and (167). Such a form of 

nonlinearity also appears in soliton dynamics and is called the 

13 

modulational nonlinearity. It is non-local (in time) as m 
<165), or is local in the form 


E = - c Z yjb- /y 1 = -Lv: N llu, - < 168 ) 

N ~ 


as obtained bs' assuming that the temporal derivative is small 
as compared to the spatial derivative. This is the case with 
the strong turbulence in cascade process. 

By taking the divergence, we transform (168) into the 
following Poisson equation 

V- E = V- -!r V ■ N u.u 

= V.M. V.ti. = A ( 169 ) 

J L L J 

The approximation is made bv assuming that the transverse modes 

1 9 

of turbulence are dominant m wave scattering. The Poisson 
equation clearly indicates the locally nonlinear character of 
the source in driving the turbulence. 

An equation of state is obtained, first in the barotropic 

A A 

form (7> for relating E to N through the speed of propagation 

A a 

(2), and then m the Poisson form <169> for relating E to u by 
assuming a transverse scattering of waves. 


l 
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9 , 2 Scattering function r* 3 

The scattering function of the macro-group is 

. « . ° -A 

A. = A n. 

° o / a o v <17 

= V i^c v < u r -(?■“'■ V ) ’ 

</ L 1 J V 6 *- <f x 

by (169), The intensity is obtained by the mean square value 



>][%•£ 


( 171 ) 


The differential 7 does not extend beyong the immediate function. 
For 




0 o 

1 ^ 





< 172 ) 


we shall decouple into double products by approximation. The 
approximation is legitimate when it does not deal with a 
process, as is the case with (170). We find 


or 


<**>**<?«< %<£> , 

<yi> ^ 2 ////^divr^r- ikUk^ik^ 

* <.#*y «;&?> (u°tk~) icrjy- 


< 173) 


< 174 ) 


by Fourier transformation, where X is the factor of truncation 
of the Fourier decomposition. 

A homogenous turbulence has the following property 




( 175a ) 


and similarly, we have 


i 
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falser))-- ftfafU <175b) 


By the properties <17 >, we can simplify < 173 > into 




4 < 1 / 6 > 

where R~ is the vorticity function as defined by < 154 >. 

By a Fourier transformation of < 168 > and with the aid of 
(ltO"), we obtain 

Ol , Col 2 

R = 4 -R ^ < 177 > 

M E ^ 

where R is the vorticity function of u -fluctuations, as 
defined by < 154 ), and 

S°1 fk 

s 2 / olk n r (* / 

( 178 ) 


if 7 s 2 r V k i km 

E Vrt 


is the vorticity function of E -fluctuations 



< 179 > 

The inner integral accentuates large values of k' near k", and 


as a first approximation can be written as 

"i 


k’Fjk")- 


( 180 ) 


yielding the relation 


fa" k- l FUl s f fa ^[Fjkvf , 


(. 1 S 1 > 


or, after a differentiation. 


3 
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F(k) • ( 182) 

The approximation is valid for small k and even smaller k” as 
in the case with the coupling subrange. 

The relations (177) and (182) are the equation of state in 
the s-space and k-space . respectively. 


10. Spectrum of the enhanced turbulence. 

10.1 Enercn ? s pectrum i n the coupling Subranne 

In the coupling subrange, the turbulent energy is enhanced 
at the espense of the potential energy. The governing equation 
of energy balance is (130), rewritten as 


or 


3 = K W K™ 

X 07 


< lS3a) 


,j\k' F £ (k')\\k') = * ( 


183b) 


C°1 


by (130). (122) and (155), where R has the dimension of 


( time) 


-2 


The equation ( 183a.) describes the balance between the 
coupling function as a spectral source and the cascade transfer 
as a spectral sink at the larger end of the spectrum. 
Therefore, the formula (153a) of the transfer by direct cascade 
is taken. 


An inspection of (183b) reveals that any function having 
the dimension of time is independent of k, 

For describing the energy flow, we differentiate (183a) 
with respect to k, as denoted by an upper dot, to get 


J 
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3 = K f '’ K M + K t] 

S K M t w 


C'] 


; 


< 184) 


in view of small R . In terms of the spectral distributions, 

<184) is, by ( 183b > , 

* f fdk" rjk"j [xfaj] 3 , 


or after a rearrangement^ 


r« r * ai r * 


-2 


( 185) 


The factor 


f t( k ) '&) 

-i 


having the dimensions of (time) , is independent of k. Another 


differentiation of (185) leads to 

-2 / -3 


W) ’ V h ■ 


The factor 


, -2 PJ / C'7, 

i •’Ti/Vj, ) 


( 186) 


( 187) 


is again independent of k. 

-C°J 


M 


Since ** is the time of steepening of sawtooth, and T is 

C At 

the modulation time for the cascade of energy, the balance 

between the two processes gives the duration t for the 

u. 

enhancement of turbulence. 

By the definitions (1.23) and (157), we can write 


J At T Z &(*)] f \ °tV Mv)j 


< 188) 


as related to the normalized second moment of the orbit 
function 




V = 0 


i 
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10.2 Spectrum of field and density fluctuations in the poudI i ng 
, ubrange 

By the use of the equation of state ( 182) , the relation 
(165b) and the spectral result (186), we derive the spectral 
distributions of field and density fluctuations as follows: 


W - 

F N Ck) -- f 


( 189a) 


( 189b) 


10.3 Spectral intensities 

By introducing the spectral intensities 



/V F(k') 


( 190a) 


( 190b> 

(190c) 


the results <1S6), ( 139a) and < 189b) for the spectral densities 

can be converted into the following : 



< 192a) 


( 192b) 


< 192C) 


i 
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It is to be noted that the spectral results for 
velocity and field fluctuations are independent of c, while the 
spectral results for density fluctuations will depend on c. For 
the density spectrum, c enters as a constant parameter, 
implying an isothermal propagation. 

The coupling subrange is located between the production 
subrange and the inertia subrange, with the mean density 
gradient and the rate of energy dissipation as the respective 
external parameters. It possesses the two parameters: the 
constant speed of propagation in isothermal gas that enters 
into the equation of propagation as en external parameter, and 
the constant time scale of enhancement of turbulence as an 
internal parameter. 

11. Discussions and conclusions 

In incompressible turbulence, the agent of production of 
turbulence is usually the mean gradient of the velocity itself. 
However, if the density may fluctuate. such as in convective 
turbulence, acoustic turbulence, density -induced turbulence, 
and other types of compressible turbulence, the turbulent 
motion may interact with the large- amplitude density waves. If 
the turbulence motion is intensified by the density waves, we 
call it density-induced turbulence. The processes are as follows: 
often a mean density gradient is more likely to be found than a 
mean velocity gradient. It excites all the nonlinear mechanisms 
that controls the density waves, i.e. the density-velocity 
coupling among large eddies. and the cascade transfer from 


j 
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large toward small eddies . 

If the coupling by large eddies is more intense than the 
transfer, the density fluctuations will mainly go to feed the 
kinetic energy for its continuing cascade toward smaller ed- 
dies, and consequently the density spectrum will fall rapidly. 
At increasing wavenumbers. the tail portion of the velocity 
spectrum is governed by the balance between the cascade trans- 
fer and the viscous dissipation. 

We need a mathematical model to describe these processes 
i n density -induced turbulence. In Section 2 s we have named 
several mathematically equivalent models. The Navier-Stokes 
model (a> and its kinetic correspondent (d) are too formal for 
the purpose. The Riemann system (c) has the merit of illustra- 
ting the build-up of the kinetic energy of turbulence at the 
expense of the wave energy. By the symmetry, the system does 
not distinguish the different physical roles from the two 
equations. The acoustic turbulence model (b) that has clearly 
put the propagation into evidence is more suitable to treat the 
density -induced turbulence. In order to include those transport 
properties which take their origin from the velocity distribu- 
tion of eddies, we choose a kinetic approach in the wave- 
kinetic model (e - *, upon replacing the Navier-Stokes equation of 
motion by a master equation. By group-scaling, we decompose the 
master equation into scaled equations for the distribution 
functions of many scales. The closure is found by considering 
those groups which form a relaxational cluster as an effective 
medium having self-consistent ly the same transport property as 
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the cluster itself. As a result, the kinetic equation for the 
macro-distribution can be used to derive the transport func- 
tions [coupling function and transfer function^*, the transport 
Properties (eddy diffusivity and eddy viscosity). and finally 
to find the spectral functions. The spectral intensities in the 
coupling subrange are found to have the power law k~— for the 
velocity and field fluctuations and the power law for the 
densitv fluctuations. This is the first time that a kinetic 
equation for the singlet distribution (macro-distribution) can 
determine the spectral structure of turbulence without the need 
of the pair-distribution. 
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